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ABSTRACT A full-atomic molecular model of human apurinic/apyrimidinic endonuclease APE1, an important 
enzyme in the DNA repair system, has been constructed. The research consisted of hybrid quantum mechanics/ 
molecular mechanics modeling of the enzyme-substrate interactions, as well as calculations of the ionization 
states of the amino acid residues of the active site of the enzyme. The choice of the APE1 mechanism with an 
Asp210 residue as a proton acceptor was validated by means of a generalization of modeling and experimental 
data. Interactions were revealed in the active site that are of greatest significance for binding the substrate and 
potential APE1 inhibitors (potential co-drugs of interest in the chemo- and radiotherapy of oncological dis- 
eases). 
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ABBREVIATIONS BER — base excision repair; AP - apurinic/apyrimidinic; APE1 — human apurinic/apyrimidinic 
endonuclease 1; QM/MM - quantum mechanics/molecular mechanics; MD - molecular dynamics. 



INTRODUCTION 

DNA damages occur frequently as a result of rep- 
lication errors or on exposure to various exo- and 
endogenic factors, such as ultraviolet radiation and 
oxidative stress. In order to ensure genomic stability, 
mammals possess enzymatic repair systems - direct 
repair, base/nucleotide excision repair, and recombi- 
nation mechanisms - which facilitate the elimination 
of most forms of DNA damage [1 — 3]. The pharmaco- 
logical inhibition of the repair systems is a promis- 
ing method for improving the efficacy of oncological 
therapy. The reason which accounts for this fact is 
that the repair systems resist the effect of the chemo- 
therapeutic agents (e.g., temozolomide and cisplatin 
[4]) which damage the DNA in order to kill the tu- 
mor cell. Therefore, the selective inhibition of the 
enzymes participating in the DNA repair processes 
can be used as an accompanying therapy. The agents 
whose binding in the enzyme active site affects the 
residues directly participating in the catalytic mecha- 
nism are reasonably expected to exhibit the highest 
levels of efficacy. Therefore, adequate data regard- 
ing the organization of the active site of the target 



enzyme, the charge distribution, and analysis of the 
interactions determining the strength of substrate- 
inhibitor binding are fundamentally necessary in the 
search for novel potential pharmaceuticals in the field 
of cancer therapy. 

Apurinic/apyrimidinic endonuclease 1 (APE1) is the 
key enzyme in the DNA repair system, known as "base 
excision repair" (BER). Apurinic/apyrimidinic (AP) 
sites are the deoxyribose residues in the DNA molecule 
without a nitrogenous base; they result from the enzy- 
matic hydrolysis of the N-glycoside bond of a damaged 
nucleotide and actually are the intermediates of the 
BER process. Furthermore, the AP sites may spontane- 
ously emerge in cells due to apurinization [5]. According 
to current estimations, up to 10,000 AP sites are formed 
in mammalian cells per day [6]. The APE1 endonuclease 
recognizes the AP sites and hydrolyzes their 5'-phos- 
phodiester bond for subsequent replacement with an 
undamaged nucleotide [5, 7]. The data obtained in labo- 
ratory and clinical studies attest to the significant role 
played by this enzyme in the development of a tumor 
and in the appearance of tumour resistance to antitu- 
mor agents [8]. 
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There exists a wide variety of viewpoints concerning 
the catalytic mechanism underlying the action of endo- 
nuclease APE1. The first crystal structure of the enzyme 
was obtained in 1997 (PDB ID lbix) [9]. In their descrip- 
tion, the authors proposed a mechanism in which the 
role of the general base in the catalysis was attributed to 
the His309 residue. In this hypothetical mechanism, the 
uncharged residue His309, jointly with Asp283, forms a 
proton relay system similar to the one formed by serine 
proteases. The only difference is that a water molecule 
acts as an activated nucleophilic agent (Fig. 1A). The role 
of a metal ion in this mechanism consists in the binding 
and polarization of the negatively charged phosphate 
group of the substrate and in the stabilization of the 
transition state of the enzymatic reaction. 

The fundamental significance of another residue in 
the active site (Asp210) for the catalysis was demon- 
strated in studies performed using site-directed muta- 
genesis: mutant forms of the enzyme with the substi- 
tutions Asp210Ala and Asp210Asn almost completely 
lost their catalytic properties (more than 25,000-fold 
reduction in activity was observed) as compared to the 
wild-type enzyme [10]. The determination of the crys- 
tal structures of human APE1 in complex with DNA 
derivatives resulted in a major revision of the assump- 
tions regarding the mechanism of action of the enzyme 
[11]. One of the ascertained structures (PDB ID lde8) 
is a complex of an inactive enzyme containing no metal 
ions with a substrate analogue, whereas the second 
structure (PDB ID lde9) contains a metal (bivalent 
manganese) ion and the enzyme-bound DNA analogue 
of the substrate after catalytic cleavage. The concep- 
tion was made regarding the structure of the enzyme- 
substrate complex, which simultaneously contains both 
the substrate analogue and a metal ion via the combina- 



tion (spatial superposition) of the structures. Although 
the resulting model structure of the enzyme— substrate 
complex does not contain water molecules potentially 
capable of attacking the substrate, the arrangement 
of the residues in the active site before and after the 
catalytic process allowed making assumptions regard- 
ing the alternative mechanism of the catalytic reaction 
[11]. In the scheme proposed, the Asp210 residue acts as 
a general base activating the water molecule, whereas 
the His309 residue, along with the metal ion, partici- 
pates in the binding and coordination to the phosphate 
group of the substrate (Fig. IB). It is assumed that the 
positive charge of the His309 residue participates in 
the catalytic process, which is presumably facilitated 
by the proximate location of the Asp238 residue. The 
authors [11] attribute the major stabilizing function in 
the formation of the transition state of the enzymatic 
reaction to the Asn212 residue. 

Hypotheses postulating that a secondary metal bind- 
ing area exists in the enzyme active site have been put 
forward in subsequent crystallographic [12] and mo- 
lecular dynamics (MD) [13] studies. The "two metal 
ions" mechanism of action of APE1 [12] (similar to that 
revealed in a related enzyme, endonuclease Endo IV) 
and the "moving metal" mechanism [13], involving 
the moving of the magnesium ion between two bind- 
ing sites during the catalytic process, were proposed in 
these works. 

It should be noted however that the NMR study uti- 
lizing the 25 Mg isotope [14] did not confirm the hypoth- 
esis of secondary magnesium ion binding in the active 
site of endonuclease APE1, thereby casting doubt on 
the "two metal ions" and the "moving metal" mecha- 
nisms. The authors of study [14] attribute the results of 
crystallographic studies [12] to the artefacts caused by 





Fig. 1 . Basic concepts of 
APE1 catalytic mecha- 
nism found in literature: 

(A) the His309 residue 
activates a water mol- 
ecule acting as a base; 

(B) the Asp2 10 residue 
activates a water mol- 
ecule acting as a base, 
while the His309 residue 
participates in substrate 
binding. The link atoms 
used in the present work 
on QM/MM modeling of 
APE1 are depicted as • . 



Asp210 



Asp210 
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the use of the lead ion instead of the magnesium ion. In 
turn, the effect of the motion of the metal ion during 
the MD modeling can be caused by the inaccuracy and 
approximations of the classical MD method. 

Mundle et al. [15] proposed a two-step variant of the 
APE1 mechanism involving the Tyrl71 residue acting 
in the form of the phenolate ion for a direct nucleophilic 
attack on the phosphate group of the substrate. This 
conclusion was drawn on the basis of the data obtained 
via site-directed mutagenesis on the 171 st position. The 
kinetic studies of the catalytic properties of the mutat- 
ed forms Tyrl71Ala, Tyrl71Phe, and Tyrl71His dem- 
onstrated a fall in enzymatic activity by almost five or- 
ders of magnitude. It should be noted however that the 
authors [16] subsequently admitted the inconsistency 
of the previously proposed two-step scheme and sup- 
ported the one-step mechanism, in which the His309 
residue acts as the general base activating the water 
molecule, while the Tyrl71 residue participates in the 
binding and proper orientation of the substrate. 

The molecular modeling methods could assist con- 
siderably in the study of the mechanism of action of 
the enzyme; however, this approach has not been ex- 
tensively explored in the study of APE1. Modeling of 
the inhibitor-enzyme interaction was performed with 
no allowance for the ionization state of the inhibitor, 
thereby making the interpretation of the obtained re- 
sults more difficult [17]. 

Thus, there exists no unambiguous conception re- 
garding the mechanism of action of human endonu- 
clease APE1 or the role of the amino acids of the active 
site in binding and catalysis. Therefore, shedding more 
light on the structure of the active site, the nature of 
the interactions between the enzyme and substrate or 
inhibitors and the involvement of the active site resi- 
dues into the catalytic mechanism of APE1 appears to 
be a topical task in molecular modeling. 

MATERIALS AND METHODS 

Molecular modeling software 

The ionization states of the amino acid residues were 
calculated using the PROPKA 2.0 software [18, 19]. 
The preparation of the initial structure for simulations 
and the trajectory analysis were performed using the 
AmberTools 1.2 package (http://ambermd.org). The 
energy minimization and MD simulations were per- 
formed using the Amber 10 package [20, 21]; the mo- 
lecular docking was performed using the Lead Finder 
1.1.14 program (MolTech Ltd., Russia) [22]. Modeling 
of the spatial structure of 6-hydroxy-DOPA was per- 
formed using the ACD/ChemSketch 8.17 program [23]. 
The visualization of the structures and trajectories was 
carried out using the VMD 1.8.6 software [24]. Paral- 



lel computations of the molecular dynamic trajectories 
were run on the SKIF-MSU "Chebyshev" supercom- 
puter (MSU Research Computing Center). 

Structure preparation 

The initial model of the APE1 enzyme-substrate com- 
plex was built on the basis of the lde8 crystallographic 
structure [11]. The coordinates of the attacking water 
molecule were calculated by docking. The coordinates 
of the manganese ion were transferred from the lde9 
structure; the metal type was subsequently substituted 
for magnesium. The structure of the enzyme-substrate 
complex was protonated and placed into a box of TIP3P 
type water molecules with the shortest distance of 
12 A between the box edge and protein. Sodium ions 
were added to neutralize the charge of the system. The 
ff99SB force field [25] was used to describe the protein 
and DNA molecules; parameters from the R.E.DD.B da- 
tabase (http:/ /q4md-forcefieldtools.org) [26] were used 
to describe the AP site. 

Energy minimization and molecular dynamics 

The model of the APE1 enzyme— substrate complex 
was equilibrated and subjected to calculation of the 
1000 ps MD trajectory according to the following pro- 
tocol. Initially, the two-stage energy minimization of 
the solvated system was performed. At the first stage 
(2,500 steps of the steepest descent algorithm followed 
by 2,500 steps of the conjugate gradient algorithm), 
the molecular mechanics description of the system was 
performed with the coordinates of the protein, DNA, 
and magnesium ion being fixed by position restraints 
fc(Ax) 2 with a force constant of 2 kcal/(mol-A 2 ). At the 
second stage (5,000 steepest-descent steps followed by 
5,000 conjugate-gradient steps), the system was divid- 
ed into a quantum mechanical (QM) region and a mo- 
lecular mechanical (MM) region; energy minimization 
was performed without any restraints. The QM region 
included the side chain of the Asp210 residue of the 
active site, the attacking water molecule, and the AP 
site fragment; this region was described by the RM1 
semi-empirical Hamiltonian [27]. The link atom model 
was used to make allowance for the bonds crossing the 
QM-MM boundaries. 

After energy minimization, using the aforemen- 
tioned QM/MM division the system was heated from 

0 to 300 K over 50 ps (with positional restraints of 

1 kcal/(molA 2 ) on the protein, DNA, and magnesium 
ion), equilibrated over 500 ps at 300 K, and finally sim- 
ulated for 1,000 ps. All simulations were performed us- 
ing periodic boundaries and the PME (Particle Mesh 
Ewald) method to allow for long-range electrostatic in- 
teractions. The cut-off radius of the non-bonded inter- 
actions was 10 A. The system was heated at a constant 
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volume; the equilibration and 1,000 ps trajectory simu- 
lation were performed under constant pressure. The 
temperature was controlled by the Langevin method. 
The integration time step was 0.002 ps. Interatomic 
distances and angles in the active site of APE1 were 
estimated by analyzing the 1,000 ps trajectory of equi- 
librium simulation. 

Molecular docking 

The model for performing the molecular docking pro- 
cedure was obtained as follows. Water molecules, so- 
dium ions, and the DNA substrate analogue molecule 
were removed from the structure of the solvated APE1 
enzyme-substrate complex after energy minimization. 
The energy grid map surrounding the AP-site area was 
then calculated. Finally, a potential inhibitor molecule, 
6-hydroxy-DOPA, was docked into the active site with 
the use of the genetic search algorithm implemented in 
the docking program. 

RESULTS AND DISCUSSION 

Ionization states of the active site residues 

Based on the results of the calculation of the ionization 
states of the active site residues using the PROPKA 2.0 
method, it was determined that the His309 residue is 
protonated (the calculated pK a value of 8.6 matches the 
pK 2 value of the experimentally determined pH pro- 
file of enzyme activity) under optimal conditions of the 
hydrolysis of the phosphodiester bond (pH 7—8 [12]), 
whereas the Asp210 residue is deprotonated (the calcu- 
lated pK a value of 6.2 is close to the pKl value of the pH 
profile of enzyme activity). Thus, it can be said that the 
deprotonated and negatively charged Asp210 residue 
acts as a general base in catalysis, whereas the posi- 
tively charged His309 residue participates in the bind- 
ing of the negatively charged phosphate group of the 
substrate and in stabilization of the reaction intermedi- 
ate product. Therefore, when building the full-atomic 
model of APE1, Asp210 and His309 were modelled in 
their charged forms. 

Model of the enzyme-substrate complex 
and the deduced catalytic mechanism 
of action of APE1 endonuclease 

The starting solvated model of the enzyme-substrate 
complex APE1 was created on the basis of the lde8 and 
lde9 crystallographic structures as described in Mate- 
rials and Methods. It was then necessary to optimize 
the atom positions within the model (especially the 
coordinates of the hydrogen atoms added); therefore, 
a two-stage minimization of the energy of the system 
was performed. Molecular-mechanical minimization 
was performed at the first stage to remove the larg- 



est strains in the system. Refinement of the active site 
structure was carried out at the second stage using the 
hybrid QM/MM method for energy minimization us- 
ing the RM1 Hamiltonian. The stability of the result- 
ing structure was confirmed by calculation of the 1,000 
ps QM/MM MD trajectory. The calculated interatomic 
distances in the APE1 active site are listed in Table. 
It was demonstrated via the analysis of the resulting 
model that substrate binding in the active site of APE1 
is accompanied by the formation of a number of bonds 
and interactions of a different nature. Among these, 
the hydrophobic interactions of deoxyribose of the AP 
site in the hydrophobic pocket formed by the Leu282, 
Phe266, and Trp280 residues should be noted. The free 
hydroxyl group in deoxyribose of the AP site also forms 
a hydrogen bond with the backbone carbonyl group 
of the Ala230 residue. The phosphate group at the 3' 
terminus of the AP site is held by the positive charge of 
the Argl77 residue. The phosphate group under attack 
electrostatically interacts with the magnesium ion and 
forms hydrogen bonds with the side chains of Asnl74, 
Asn212, and His309. The hydroxyl group of the Tyrl71 
residue is oriented towards the oxygen atom of the 
leaving group. 

The orientation of the attacking water molecule 
is provided via the interaction with the general base 
Asp210, carbonyl group of Asn212, and the side chain of 
His309; however, the interactions H 2 0:H2 - Asn212:ODl 
and H 2 0:0 - His309:HE2 are not conventional hydrogen 
bonds, since the average 0-H2-OD1 and NE2-HE2-0 
angles are 137° and 122°, respectively. In the case of a hy- 
drogen bond, these values should to be at least 150°. The 
reactive conformation of the carboxyl group of the gen- 
Table. Distance characteristics of the APE1 enzyme- 
substrate complex obtained via the equilibrium 1 ,000 ps 
QM/MM MD simulation. Mean values are presented 
together with the standard deviation 



Interaction 


Distance, A 


H 2 0:0-AP site:P 


1.91 ± 0.03 


H 2 O:O-His309:HE2 


2.52 ± 0.17 


H 2 O:H 1 -Asp210:ODl 


1.49 ± 0.07 


H 2 0:H 2 -Asn212:ODl 


2.62 ± 0.47 


AP site:01P-Mg 2+ 


1.84 ± 0.04 


AP site:01P-His309:HE2 


1.78 ±0.09 


AP site:02P-Asn212:HD2 


2.08 ± 0.32 


dC5:03'-Mg 2+ 


1.95 ±0.06 


Asp210:OD2-Asn212:H 


2.22 ± 0.20 
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Fig. 2. The active site of the full-atomic model of APE1 
enzyme-substrate complex. Hydrogen bonds are shown 
as dotted lines. Electrostatic interactions other than 
hydrogen bonds are shown as dashed lines. The arrow 
shows the direction of nucleophilic attack. 

eral base Asp210 in the enzyme is maintained via the in- 
teraction between its side chain and the backbone amino 
group of Asn212 (Fig. 2). 

Oriented and polarized by Asp210, His309, and the 
metal ion, the water molecule can attack the phosphate 
group of the substrate with the simultaneous trans- 
mission of the proton to the general base, the Asp210 
residue. The intermediate resulting from the attack is 
shaped as a trigonal bipyramid and stabilized via the 
following interactions in the enzyme active site: oxy- 
gen atoms in the "apexes" of the bipyramid interact 
with the side chains of the His309 and Tyrl71 residues; 
the trigonal base of the bipyramid is placed between 
Asnl74, Asn212, and the magnesium ion {Fig. 3). 

The data pertaining to the mutagenesis on Tyrl71 [15] 
attest to the crucial role of this residue in the catalytic 
mechanism of APE1; however, in contrast to the previ- 
ous suggestions [15, 16], we assume a different role for 
this residue in the catalysis. The proximate location of 
the positive charge of magnesium ion and the Argl56 
residue are supposed to facilitate proton migration from 
the hydroxyl group of Tyrl71. Therefore, we consider 
this residue to be a potential proton donor for the leav- 
ing group, which is a strong base. The weaker influence 
of the mutations at position 171 on substrate binding in 
comparison with the decrease in the catalytic constant 
[15] is attributed to the weak interaction between the 



residue and the substrate at the earlier stages of the 
reaction preceding the catalytic process, which fully 
matches the assumption made. In the catalytic trans- 
formation, the leaving group apparently approaches the 
chain of the Tyrl71 residue, enabling proton transfer. 

As the reaction proceeds, the less stabilized P-O" 
bond located at the base of the bipyramid and directed 
towards the Asn212 residue is subsequently trans- 
formed into a double P=0 bond. The P— O bond di- 
rected towards the Tyrl71 residue is simultaneously 
broken; as a result, the leaving group takes the proton 
away from the hydroxyl group of tyrosine {Fig. 3). 

Restoration of the catalytically active site (deproto- 
nation of the general base Asp210 and protonation of 
the acid Tyrl71) occurs via the interaction with sur- 
rounding water molecules. 

Analysis of the binding of substrate and 
mechanism-dependent inhibitors 

A large number of charges and polar groups participate 
in the substrate binding and intermediate stabilization 
during the reaction. The character of these interactions 
and the ionization states of the amino acid residues of 
the active site define certain requirements to the struc- 
ture of the substances capable of binding to the active 
site of APE1. It is necessary to implement the most sig- 
nificant interactions when constructing efficient inhibi- 
tors of the enzyme. The existence of the hydrophobic 




Fig. 3. The schematic picture of the hypothetical structure 
of the transition state (intermediate shaped as a trigonal 
bipyramide) and its conversion upon a hydrolysis reaction 
catalyzed by endonuclease APE1 . 
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Fig. 4. Positions 
of the substrate 

(A) and the 
potential inhibitor 

(B) in the ac- 
tive site of the 
full-atomic APE1 
model. Residues 
of the hydropho- 
bic pocket are 
shown in yellow: 
Phe266, Trp280, 
Leu282. 



binding site, along with the number of polar groups and 
groups with various charges, further complicates the 
search for low-molecular-weight compounds of appro- 
priate structure. Amino acids are a class of natural com- 
pounds the structure of which simultaneously contains 
substituents of a different nature, capable of both per- 
forming the hydrophobic and electrostatic interactions 
and of acting as hydrogen bond donor or acceptor. It 
has been noted [17] that 6-hydroxy-DOPA is a potential 
inhibitor of APElendonuclease activity; however, an 
incorrect allowance was made for the ionization state 
of the inhibitor when modeling the enzyme-inhibitor 
interaction, making it more difficult to interpret the re- 
sults. To clarify enzyme interactions with the potential 
inhibitors of this structure, molecular docking of dif- 
ferent amino acids and their derivatives into the APE1 
active site was carried out, taking into account the ion- 
ization states of the inhibitors and active site residues. 

It was demonstrated via the analysis of the molecu- 
lar docking results that the presence of the carboxyl 
group allows the selected compounds to bind to the 
metal ion and the His309 residue, whereas the hydro- 
phobic substituent (e.g., phenyl radical) can occupy the 
hydrophobic pocket {Fig. 4). The introduction of hy- 
droxyl substituents to the phenyl radical may lead to 
the formation of additional hydrogen bonds with the 
polar residues of the enzyme active site. Furthermore, 
the binding of the inhibitor to the charged residue of 
the general base is one of the factors presumably de- 
termining the inhibition efficiency. 

CONCLUSIONS 

The aim of this work was to select the most reliable 
mechanism of action of APE1 on the basis of a mo- 
lecular modeling, analysis of the full-atomic model of 
APE1, and a critique of the experimental results and 



assumptions previously made in the literature. For this 
purpose, calculation of the ionization states of the ac- 
tive site residues and a hybrid QM/MM modeling of 
the enzyme-substrate complex containing a water 
molecule capable of attacking the substrate were car- 
ried out. As a result of the investigation conducted, it 
was demonstrated that the Asp210 is likely to act as 
the general base in the catalytic mechanism, whereas 
the His309 residue, being protonated (and positively 
charged), participates in the binding of the phosphate 
group of the substrate. The analysis of the molecular 
dynamic trajectory of the enzyme-substrate complex 
attested to its high reactivity and confirmed the valid- 
ity of the molecular modeling performed. 

The most important interactions in the active site de- 
termining the efficiency of binding of the substrate and 
the potential enzyme inhibitors (which are promising 
co-drugs of interest in the chemo- and radiotherapy 
of oncological diseases) were revealed. An assumption 
regarding the role of the Tyrl71 of the active site of 
APE1 as the residue capable of ceding the proton to 
the leaving group of the substrate was made. Thus, the 
investigation enabled to establish a consistent mecha- 
nism of action of the enzyme. Furthermore, it allowed 
to summarize MD data, as well as the experimental re- 
sults of kinetic studies and the other published data. 
In the next step, we plan to use higher level QM/MM 
methods to calculate the energy barrier of the reaction 
catalyzed by endonuclease APE1, in compliance with 
the mechanism proposed, and to screen for effective 
inhibitors with the use of the constructed mechanistic 
full-atomic model of the enzyme. • 

This work was supported by the Ministry of Education 
and Science of the Russian Federation ( Government 
contract M 16.512.11.2240). 



VOL. 4 JVs2(13) 2012 | ACTA NATURAE | 85 



RESEARCH ARTICLES 



REFERENCES 

1 . Lindahl T, Wood R.D. // Science. 1999. V. 286. P. 1897-1905. 

2 . Hoeijmakers J.H. // Nature. 2011. V. 411. P. 366-374. 

3 . Rechkunova N.I., Krasikova Y.S., Lavrik O.I. // Biochem- 
istry (Moscow). 2011. V. 76. P. 24 - 35. 

4 . Abbotts R., Madhusudan S. // Cancer Treatment Rev. 
2010. V. 36. P. 425-435. 

5 . Nevinsky G.A. // Biochemistry (Moscow). 2011. V. 76. P. 
94 - 117. 

6 . Wilson D.M., Simeonov A. // CMLS. 2010. V. 67. P. 3621- 
3631. 

7 . Dyrkheeva N.S., Khodyreva S.N, Lavrik O.I. // Molecular 
biology (Moscow). 2007. V. 41. P. 402 - 416. 

8 . Fishel M.L., Kelley M.R. // Mol. Aspects Medicine. 2007. V. 
28. P. 375-395. 

9 . Gorman M.A., Morera S., Rothwell D.G., de La Fortelle E., 
Mol CD., Tainer J.A., Hickson I.D., Freemont PS. // EMBO 
J. 1997. V. 16. P. 6548-6558. 

10 . Erzberger J.P, Wilson D.M. // J. Mol. Biol. 1999. V. 290. P. 
447-457. 

11 . Mol CD., Izumi T, Mitra S., Tainer J.A. // Nature. 2000. 
V. 403. P. 451-456. 

12 . Beernink P.T., Segelke BW, Hadi M.Z., Erzberger J.P, 
Wilson D.M., Rupp B. // J. Mol. Biol. 2001. V. 307. P. 1023- 
1034. 

13 . Oezguen N., Schein C.H., Peddi S.R., Power T.D., Izumi 
T, Braun W. // Proteins: Structure, Function and Bioinfor- 
matics. 2007. V. 68. P. 313-323. 

14 . Lipton A.S., Heck RW, Primak S., McNeill D.R., Wilson 
D.M., Ellis PD // J. Amer. Chem. Soc. 2008. V. 130. P. 9332- 
9341. 



15 . Mundle ST., Fattal M.H., Melo L.F, Coriolan J.D., 
O'Regan N.E., Strauss PR. // DNA Repair. 2004. V. 3. P. 
1447-1455. 

16 . Mundle ST., Delaney J.C., Essigmann J.M., Strauss PR. 
// Biochemistry. 2009. V. 48. P. 19-26. 

17 . Simeonov A., Kulkarni A., Dorjsuren D., Jadhav A., Shen 
M., McNeill DR., Austin CP, Wilson D.M. Ill // PloS One. 
2009. V. 4. P. e5740. 

18 . Li H., Robertson A.D., Jensen J.H. // Proteins. 2005. V. 61. 
P. 704-721. 

19 . Bas DC, Rogers DM., Jensen J.H. // Proteins. 2008. V. 73. 
P. 765-783. 

20 . Case D.A., Darden T.A., Cheatham T.E. Ill, Simmerling 
C.L., Wang J., Duke R.E., Luo R., Crowley M., Walker R.C., 
Zhang W, et al. // AMBER 10. University of California. San 
Francisco. 2008. 

21 . Walker R.C, Crowley M.F, Case DA. // J. Comput. Chem. 
2008. V. 29. P. 1019-1031. 

22 . Stroganov O.V., Novikov F.N., Stroylov V.S., Kulkov V., 
Chilov GG // J. Chem. Inf. Model. 2008. V. 48. P. 2371-2385. 

23 . ACD/ChemSketch Freeware, version 8.17. Advanced 
Chemistry Development, Inc., Toronto, ON, Canada, www. 
acdlabs.com. 2005. 

24 . Humphrey W, Dalke A., Schulten K. // J. Mol. Graphics. 
1996. V. 14. M> 1. P. 33-38. 

25 . Hornak V., Abel R., Okur A., Strockbine B., Roitberg A., 
Simmerling C. // Proteins. 2006. V. 65. P. 712-725. 

26 . Chen J., Dupradeau F.-Y., Case D.A., Turner C. J., Stubbe 
J. // Biochemistry. 2007. V. 46. P. 3096-3107. 

27 . Rocha G.B., Freire R.O., Simas A.M., Stewart J. J.P. // J. 
Comp. Chem. 2006. V. 27. P. 1101-1111. 



86 | ACTA NATURAE | VOL.4 JY°2(13) 2012 



